Constraining the High-Density Behavior of Nuclear Symmetry Energy with the Tidal 

Polarizability of Neutron Stars 
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Using a set of model equations of state satisfying the latest constraints from both terrestrial 
nuclear experiments and astrophysical observations as well as state-of-the-art nuclear many-body 
calculations of the pure neutron matter equation of state, the tidal polarizability of canonical neu- 
tron stars in coalescing binaries is found to be a very sensitive probe of the high-density behavior 
of nuclear symmetry energy which is among the most uncertain properties of dense neutron-rich 
nucleonic matter. Moreover, it changes less than ±10% by varying various properties of symmetric 
nuclear matter and symmetry energy around the saturation density within their respective ranges 
of remaining uncertainty. 
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Introduction. — Understanding the nature of neutron- 
rich nucleonic matter is a major thrust of current re- 
search in both nuclear physics and astrophysics [l|. To 
reaUze this goal, many experiments and observations are 
being carried out or proposed using a wide variety of 
advanced new facilities, such as. Facilities for Rare Iso- 
tope Beams (FRIB), X-ray satellites and gravitational 
wave detectors. Most critical to interpreting results 
of these experiments and observations is the equation 
of state (EOS) of neutron-rich nucleonic matter, i.e., 
E{p,a) = Eo{p) + S{p)a^ + ©(a"), where E{p,a) and 
Eq (p) are the specific energy in asymmetric nuclear mat- 
ter of isospin asymmetry a = {pn — Pp)/p and symmetric 
nuclear matter (SNM), respectively, and S{p) is the sym- 
metry energy encoding the energy cost of converting all 
protons in SNM to neutrons. Thanks to the continuing 
efforts of both the nuclear physics and astrophysics com- 
munity over several decades, the EOS of SNM around the 
saturation density po has been well constrained. More- 
over, combining information from studying the collec- 
tive flow and kaon production in relativistic heavy-ion 
collisions in several terrestrial nuclear physics laborato- 
ries and the very recent discovery of the maximum 
mass of neutron stars Q, the EOS of SNM has been lim- 
ited in a relatively small range up to about 4.5po- The 
symmetry energy S{p) is a vital ingredient in describing 
the structure of rare isotopes and their reaction mech- 
anisms. It also determines uniquely the proton fraction 
and thus the cooling mechanism, appearance of hyperons 
and possible kaon condensation in neutron stars. More- 
over, it affects significantly the structure, such as the 
radii, moment of inertia and the core-crust transition 
density, as well as the frequencies and damping times 
of various oscillation modes of neutron stars, see, e.g., 
Refs. H, [HI for reviews. Intensive efforts devoted to con- 
straining S (p) using various approaches have recently led 
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exceptions, although the error bars for L from different 
approaches may vary broadly [6l-[l^. On the other hand, 
the high-density behavior of S{p) remains very uncer- 
tain despite its importance to understanding what hap- 
pens in the core of neutron stars and in reac- 
tions with high energy radioactive beams [isj . The pre- 
dictions for the high-density behavior of the symmetry 
energy from all varieties of nuclear models diverge dra- 
matically 13" 21 1, with some models predicting very stiff 
symmetry energies that increase continuously with den- 
sity [2l|-|26[, and others predicting relatively soft ones, or 
an S{p) that first increases with density, then saturates 
and starts decreasing with increasing density fiol 27- 
41|. These uncertainties can be traced to our poor 



knowledge about the isospin dependence of strong in- 
teraction in dense neutron-rich medium, particularly the 
spin-isospin dependence of three-body and many-body 
forces, the short-range behavior of nuclear tensor force 
and the isospin dependence of nucleon-nucleon correla- 
tions in dense medium, see, e.g. Refs. (i^. El]- Little 
experimental progress has been made in constraining the 
high density S{p) partially because of the lack of sen- 
sitive pro bes. While several observables have been pro- 
posed [l8| and some indications of the high-density S{p) 
have been reported recently [44, 45|, conclusions based on 
terrestrial nuclear experiments remain controversial (46j . 
To our best knowledge, the only astrophysical probe for 
high density S{p) proposed so far is the late time neu- 
trino signal from a core collapse supernova 47[. In this 
Letter, we show that the tidal polarizability of canonical 
neutron stars in the coalescing binaries is a very sensitive 
probe of the high-density behavior of nuclear symmetry 
energy independent of the remaining uncertainties of the 
SNM EOS and S{p) near saturation density. 
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Stellar structure and tidal polarizability. — Coalescing 
binary neutron stars are among the most promising 
sources of gravitational waves (GW). One of the most 
important features of binary mergers is the tidal deforma- 
tion neutron stars undergo as they approach each other 
prior to merger, the strength of which can give us pre- 
cious information about the neutron-star matter EOS 
[iil - IHgl l . At the early stage of an inspiral tidal effects may 
be effectively described through the tidal polarizability 
parameter A [isl. Isil - IHst defined via Qij = where 
Qij is the induced quadrupole moment of a star in binary, 
and £ij is the static external tidal field of the companion 
star. The tidal polarizability can be expressed in terms of 
the dimensionless tidal Love number ^2 and the neutron 
star radius i? as A = 2fc2-R^/(3G). The tidal Love number 



^2 is found using the following expression [49|, |5J| : 
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where Rg = 2M is the Schwarzschild radius of the star, 
and yn = y{R) can be calculated by solving the following 
first-order differential equation: 



ar 



with 



F(r) = 



r - 47rr3 {£{r) - P{r)) 
r - 2M(r) 



(2) 



(3) 



Q{r) 



47rr(5g(r)+9P(r) + /;g+gff) 



6 

47rr2 



2M(r) 



M{r) + A-nr^P{r) 
r2 (1 - 2M{r)/r) 



(4) 



The Eq. ([2) must be integrated together with the 
Tolman-Oppenhcimer-VolkofF (TOV) equation. Given 
the boundary conditions in terms of y(0) = 2, P(0) = Pc 
and M(0) = 0, the tidal Love number can be obtained 
once an EOS is supplied. Previous studies have used 
both polytropic EOSs and several popular nuclear EOSs 
available in the literature iJH^- While other parti- 



cles may be present, for the purpose of this work, it is 
sufficient to assume that neutron stars consist of only 
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FIG. 1: (Color online) The EOS of SNM and PNM as well as 
the symmetry energy as a function of density obtained within 
the lU-FSU RMF model and the SHF approach using the 
SklU-FSU parameter set. 



neutrons (n), protons (p), electrons (c) and muons (/i) in 
/3-equilibrium. 

Constrained EOS of neutron-rich nuclear matter. — 
We use two classes of nuclear EOSs within the Relativis- 
tic Mean Field (RMF) model and the Skyrme Hartre- 
Fock (SHF) approach. All of these EOSs are adjusted to 
satisfy the following four conditions within their respec- 
tive uncertain ranges: (1) reproducing the PNM EOS at 
sub-saturation densities predicted by the latest state-of- 
the-art microscopic nuclear many-body theories 2^, 60l - 
IgH I: (2) predicting correctly saturation properties of sym- 
metric nuclear matter, i.e, nucleon binding energy B = 
-16 ± 1 MeV and incompressibility Kq = 230 ± 20 MeV 
and nucleon effective mass Mq = 0.61±0.03 M at satura- 



tion density po = 0.155±0.01 fm~ ; (3) predicting a fidu- 
cial value of symmetry energy S{2pq/2>) = 26 ±0.5 MeV, 
J = S{pq) = 31 ± 2 MeV and the density slope of sym- 
metry energy L = 50 ± 10 MeV; (4) passing through the 
terrestrial constraints on the EOS of SNM between 2po 
and 4.5po and giving a maximum mass of neutron 
stars of about 2Mq assuming they are made of only the 
npep, matter without considering other degrees of free- 
dom or invoking any exotic mechanism 0, [6a | . As an ex- 
ample, two such EOSs obtained using the lU-FSU RMF 
model 1671 and the SHF using the SklU-FSU parame- 
ter set [12[ are shown in Fig. 1. By design, they both 
have the same EOS for SNM and PNM around and be- 
low pq. Thus, at sub-saturation densities the values of 
S(p) which is approximately the difference between the 
EOSs for PNM and SNM are almost identical for the two 
models. However, the values of S{p) are significantly dif- 
ferent above about 1.5po with the lU-FSU leading to a 
much stifFcr S{p) at high densities. More quantitatively, 
the S{p) with lU-FSU is 40 - 60% higher in the density 
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FIG. 2: (Color online) Energy per nucleon as a function of 
the Fermi momentum for PNM for selected models described 
in the text. 



FIG. 4: (Color online) Percentage changes in the tidal po- 
larizability of a 1.4 solar mass neutron star by individually 
varying properties of nuclear matter Ko (a), M* (b), L (c), 
and the ^ parameter (d) of the RMF model with respect to 
the value using the base model. 




FIG. 3: (Color online) The pressure of SNM given as the 
function of baryon density. Here po is the nuclear matter 
saturation density and the shaded area represents the EOS 
extracted from the analysis of 0). 

range of p/po = 3 — 4 expected to reach in the core of 
canonical neutron stars. 

To test the sensitivity of the tidal polarizability to 
variations of properties of neutron-rich nuclear matter 
around po within the constraints listed above, we build 
17 RMF parameterizations by systematically varying 
the values of Kq, Mq, L, and the ( parameter of the 
RMF model that controls the omega-meson self interac- 
tions and subsequently the high-density component 
of the EOS of SNM. Besides the constraints listed above, 
all parameter sets can correctly reproduce the experi- 
mental values for the binding energy and charge radius 
of 208 pi^ g^j^jj ^i^Q ground state properties of other closed 



shell nuclei within 2% uncertainty [69|. As a reference 
for comparisons, we select Kq = 230 MeV, Mq = 0.61 
M, i = 50 MeV, and C = 0.025 for our base model, 
which predicts po = 0.1524 fm^^^ ^ ^ -16.33 MeV 
and J = 31.64 MeV. The representative model EOSs for 
PNM at sub-saturation densities and those for SNM at 
supra-saturation densities are compared with their con- 
straints in Fig. 2 and Fig. 3, respectively. It is seen that 
the SklU-FSU and all the RMF models with 42 < L < 58 
MeV can satisfy the PNM EOS constraint. Also, they 
can all satisfy simultaneously the high density SNM EOS 
constraint with 0.02 < C < 0.03. Moreover, they all give 
a maximum mass for neutron stars between 1.94Mq and 
2.07Mq and radii between 12.33 km and 13.22 km for 
canonical neutron stars (67l | consistent with existing ob- 
servations (sl. [HgI [toI . 

Results and discussions. — First, we examine sen- 
sitivities of the tidal polarizability A of a 1.4Mq neu- 
tron star to the variations of SNM properties and the 
slope of the symmetry energy around po in Fig. 4. The 
changes of A relative to the values for our base RMF 
model are shown for the remaining RMF EOSs. It 
is very interesting to see that the tidal polarizability 
is rather insensitive to the variation of L although it 
changes up to ±10% with Kq, M* and ( within their 
individual uncertain ranges. While the averaged mass 
is M = 1.33 ± 0.05 M0, neutron stars in binaries have 
a broad mass distribution [7lj. It is thus necessary to 
investigate the mass dependence of the tidal polarizabil- 
ity. Whereas what can be measured for a neutron star 
binary of mass Mi and M2 is the mass-weighted polar- 
izability A = [Ai (Afi + I2M2) /Mi 1 O 2] /26 (El, for 
the purpose of this study it is sufficient to consider bi- 
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FIG. 5: (Color online) Tidal polarizability A of a single neu- 
tron star as a function of neutron-star mass for a range of 
EOS that allow various stiffness of symmetry energies. A 
crude estimate of uncertainties in measuring A for equal mass 
binaries at a distance of D = 100 Mpc is shown for the Ad- 
vanced LIGO (shaded light-grey area) and the Einstein Tele- 
scope (shaded dark-grey area). 



naries consisting of two neutron stars with equal masses. 
What can we learn from the tidal polarizability of light 
and massive neutron stars, respectively? Shown in Fig. 5 
are the tidal polarizability A as a function of neutron-star 
mass for a range of EOSs. Most interestingly, it is seen 
that the lU-FSU and SklU-FSU models which are dif- 
ferent only in their predictions for the nuclear symmetry 
energy above about 1.5po as shown in Fig. 1 lead to sig- 
nificantly different A values in a broad mass range from 
0.5 to 2 Mq. More quantitatively, a 41% change in A 
from 2.828 x 10^6 (lU-FSU) to 1.657 x lO^*^ (SklU-FSU) 
is observed for a canonical neutron star of 1.4 Mq. For 
a comparison, we notice that this effect is as strong as 
the symmetry energy effect on the late time neutrino flux 
from the cooling of proto-neutron stars 47| . Moreover, it 
is shown that the variation of L has a very small effect on 
the tidal polarizability A of massive neutron stars, which 
is consistent with the results shown in Fig. 4. On the 
other hand, the L parameter affects significantly the tidal 
polarizability of neutron stars with M < 1.2M0. These 
observations can be easily understood. From Eq. ([T]) the 
Love number ^2 is essentially determined by the com- 
pactness parameter M/R and the function y{R)- Both 
of them are obtained by integrating the EOS all the way 
from the core to the surface. Since the saturation den- 
sity approximately corresponds to the central density of 
a O.3M0 neutron star, one thus should expect that only 
the Love number of low-mass neutron stars to be sensi- 
tive to the EOS around the saturation density. However, 
for canonical and more massive neutron stars, the cen- 
tral density is higher than 3 — 4po, and therefore both 



the compactness M/R and y{R) show stronger sensitiv- 
ity to the variation of EOS at supra-saturation densities. 
Since all the EOSs for SNM at supra-saturation densities 
have already been constrained by the terrestrial nuclear 
physics data and required to give a maximum mass about 
2Mq for neutron stars, the strongest effect on calcula- 
tions of the tidal polarizability of massive neutron stars 
should therefore come from the high-density behavior of 
the symmetry energy. 

It has been suggested that the Advanced LIGO- Virgo 
detector may potentially measure the tidal polarizability 
of binary neutron stars with a moderate accuracy. Are 
the existing or planned GW detectors sensitive enough to 
measure the predicted effects of high-density symmetry 
energy on the tidal polarizability? To answer this ques- 
tion, as an example we estimate uncertainties in mea- 
suring A for equal mass binaries at an optimally-oriented 
distance of Z? ~ 100 Mpc 1531. 17211 using the same ap- 
proach as detailed in Rcfs. [53|, |59[. These are shown for 
the Advanced LIGO- Virgo (shaded light-grey area) and 
the Einstein Telescope (shaded dark-grey area) in Fig. 
5. It is seen that the Advanced LIGO- Virgo will un- 
likely constrain the EOS and symmetry energy at supra- 
saturation densities within the estimated uncertainty, al- 
though it is possible that a rare but nearby binary sys- 
tem may be found and provide a much more tighter con- 
straint [53|. Nevertheless, measurements for binaries con- 
sisting of light neutron stars can still help further con- 
strain the symmetry energy around the saturation den- 
sity. On the other hand, it is exciting to see that the 
narrow uncertain range for the proposed Einstein Tele- 
scope will enable it to tightly constrain the symmetry 
energy especially at high densities. 

Conclusions. — Using the EOSs for neutron-rich nucle- 
onic matter satisfying the latest constraints from both 
terrestrial nuclear experiments and astrophysical obser- 
vations, as well as the state-of-the-art nuclear many-body 
calculations for PNM EOS, we found that the tidal polar- 
izability of canonical neutron starts in coalescing binaries 
is very sensitive to the high-density behavior of nuclear 
symmetry energy, but little affected by the variations of 
SNM EOS and symmetry energy around the saturation 
density within their remaining uncertainty ranges. Fu- 
ture measurements of the tidal polarizability of neutron 
stars using the proposed Einstein Telescope will help con- 
strain stringently the high-density behavior of nuclear 
symmetry energy, and thus the nature of dense neutron- 
rich nucleonic matter. 

Acknowledgments. — We would like to thank Dr. 
J. Piekarewicz and Dr. Jun Xu for various useful 
discussions. This work is supported in part by the 
National Aeronautics and Space Administration under 
grant NNX11AC41G issued through the Science Mission 
Directorate, and the National Science Foundation under 
Grants No. PHY-0757839, No. PHY-1062613 and No. 
PHY-1068022. 



5 



* Electronic address: I Farrooh . Fat t oyev@t am uc . edu| 
^ Electronic address: Jose. Carvajal002Qmym dc.net | 

* Electronic address: William. Newton@tamuc.eduJ 
^ Electronic address: Bao-An.LiQtamuc.cdu 

[1] The 2007 US NSAC Long Range Plan, URL 
^ttp : //www. er . doe . gov/np/nsac/nsac.html 

[2] P. Danielewicz, R. Lacey, and W. G. Lynch, Science 298, 
1592 (2002). 

[3] P. B. Demorest, T. Pennucci, S. M. Ransom, M. S. E. 

Roberts, and J. W. T. Hessels, Nature 467, 1081 (2010). 
[4] J. M. Lattimer and M. Prakash, Phys. Rept. 333, 121 

(2000). 

[5] J. M. Lattimer and M. Prakash, Science 304, 536 (2004). 
[6] C. Xu, B.-A. Li, and L.-W. Chen, Phys. Rev. C82, 
054607 (2010). 

[7] W. G. Newton, M. Gearheart, and B.-A. Li (2011), 
1110.4043. 

[8] A. W. Steiner and S. Gandolfi, Phys. Rev. Lett. 108, 

081102 (2012). 
[9] J. M. Lattimer and Y. Lim (2012), 1203.4286. 
[10] M. Dutra, O. Lourengo, J. S. Sa Martins, A. Delfino, J. R. 

Stone, and P. D. Stevenson, Phys. Rev. C85, 035201 

(2012). 

[11] M. Tsang et al, Phys. Rev. C86, 015803 (2012). 

[12] F. J. Fattoyev, W. G. Newton, J. Xu, and B.-A. Li, Phys. 

Rev. C86, 025804 (2012). 
[13] M. Kutschera, Phys. Lett. B340, 1 (1994). 
[14] S. Kubis and M. Kutschera, Acta Phys. Polon. B30, 2747 

(1999). 

[15] S. Kubis and M. Kutschera, Nucl. Phys. A720, 189 
(2003). 

[16] D.-H. Wen, B.-A. Li, and L.-W. Chen, Phys. Rev. Lett. 

103, 211102 (2009). 
[17] H. K. Lee, B.-Y. Park, and M. Rho, Phys. Rev. C83, 

025206 (2011). 

[18] B.-A. Li, L.-W. Chen, and C. M. Ko, Phys. Rept. 464, 
113 (2008). 

[19] B. A. Brown, Phys. Rev. Lett. 85, 5296 (2000). 

[20] A. Szmaglinski, W. Wojcik, and M. Kutschera, Acta 

Phys. Polon. B37, 277 (2006). 
[21] A. E. L. Dieperink, Y. Dewulf, D. Van Neck, M. Waro- 

quier, and V. Rodin, Phys. Rev. C68, 064307 (2003). 
[22] A. W. Steiner, M. Prakash, J. M. Lattimer, and P. J. 

Ellis, Phys. Rept. 411, 325 (2005). 
[23] C.-H. Lee, T. T. S. Kuo, G. Q. Li, and G. E. Brown, 

Phys. Rev. C57, 3488 (1998). 
[24] C. J. Horowitz and J. Piekarewicz, Phys. Rev. Lett. 86, 

5647 (2001). 

[25] L.-W. Chen, C. M. Ko, and B.-A. Li, Phys. Rev. C76, 
054316 (2007). 

[26] Z. H. Li, U. Lombardo, H.-J. Schulze, W. Zuo, L. W. 

Chen, and H. R. Ma, Phys. Rev. C74, 047304 (2006). 
[27] V. R. Pandharipande and V. K. Garde, Phys. Lett. B39, 

608 (1972). 

[28] B. Friedman and V. R. Pandharipande, Nucl. Phys. 

A361, 502 (1981). 
[29] R. B. Wiringa, V. Fiks, and A. Fabrocini, Phys. Rev. 

C38, 1010 (1988). 
[30] N. Kaiser, S. Fritsch, and W. Weise, Nucl. Phys. A697, 

255 (2002). 

[31] P. G. Krastev and F. Sammarruca, Phys. Rev. C74, 



025808 (2006). 

[32] E. Chabanat, J. Meyer, P. Bonche, R. Schaeffer, and 
P. Haensel, Nucl. Phys. A627, 710 (1997). 

[33] J. R. Stone, J. C. Miller, R. Koncewicz, P. D. Stevenson, 
and M. R. Strayer, Phys. Rev. C68, 034324 (2003). 

[34] L.-W. Chen, C. M. Ko, and B.-A. Li, Phys. Rev. C72, 
064309 (2005). 

[35] J. Decharge and D. Gogny, Phys. Rev. C21, 1568 (1980). 
[36] C. B. Das, S. D. Gupta, C. Gale, and B.-A. Li, Phys. 

Rev. C67, 034611 (2003). 
[37] D. T. Khoa, W. von Oertzen, and A. A. Ogloblin, Nucl. 

Phys. A602, 98 (1996). 
[38] D. N. Basu, P. R. Chowdhury, and C. Samanta, Nucl. 

Phys. A811, 140 (2008). 
[39] W. D. Myers and W. J. Swiatecki, Acta Phys. Polon. 

B26, 111 (1995). 
[40] S. Banik and D. Bandyopadhyay, J. Phys. G26, 1495 

(2000). 

[41] P. R. Chowdhury, D. N. Basu, and C. Samanta, Phys. 

Rev. C80, 011305 (2009). 
[42] C. Xu and B.-A. Li, Phys. Rev. C81, 064612 (2010). 
[43] C. Xu, A. Li, and B.-A. Li (2012), 1207.1639. 
[44] Z. Xiao, B.-A. Li, L.-W. Chen, G.-C. Yong, and 

M. Zhang, Phys. Rev. Lett. 102, 062502 (2009). 
[45] P. Russotto et al., Phys. Lett. B697, 471 (2011). 
[46] W. Trautmann and H. H. Wolter, Int. J. Mod. Phys. 

E21, 1230003 (2012). 
[47] L. F. Roberts, G. Shen, V. Cirigliano, J. A. Pons, 

S. Reddy, and S. E. Woosley, Phys. Rev. Lett. 108, 

061103 (2012). 

[48] E. E. Flanagan and T. Hinderer, Phys. Rev. D77, 021502 

(2008) . 

[49] T. Hinderer, Astrophys. J. 677, 1216 (2008). 

[50] T. Binnington and E. Poisson, Phys. Rev. D80, 084018 

(2009) . 

[51] T. Damour and A. Nagar, Phys. Rev. D80, 084035 

(2009) . 

[52] T. Damour and A. Nagar, Phys. Rev. D81, 084016 

(2010) . 

[53] T. Hinderer, B. D. Lackey, R. N. Lang, and J. S. Read, 

Phys. Rev. D81, 123016 (2010). 
[54] S. Postnikov, M. Prakash, and J. M. Lattimer, Phys. Rev. 

D82, 024016 (2010). 
[55] L. Baiotti, T. Damour, B. Giacomazzo, A. Nagar, and 

L. RezzoUa, Phys. Rev. Lett. 105, 261101 (2010). 
[56] L. Baiotti, T. Damour, B. Giacomazzo, A. Nagar, and 

L. RezzoUa, Phys. Rev. D84, 024017 (2011). 
[57] B. D. Lackey, K. Kyutoku, M. Shibata, P. R. Brady, and 

J. L. Friedman, Phys. Rev. D85, 044061 (2012). 
[58] F. Pannarale, L. RezzoUa, F. Ohme, and J. S. Read, 

Phys. Rev. D84, 104017 (2011). 
[59] T. Damour, A. Nagar, and L. Villain, Phys. Rev. D85, 

123007 (2012). 

[60] A. Schwenk and C. J. Pethick, Phys. Rev. Lett. 95, 
160401 (2005). 

[61] E. N. E. van Dalen and H. Muther, Phys. Rev. C80, 
037303 (2009). 

[62] K. Hebeler and A. Schwenk, Phys. Rev. C82, 014314 
(2010). 

[63] S. Gandolfi, A. Y. lUarionov, S. Fantoni, F. Pederiva, and 
K. E. Schmidt, Phys. Rev. Lett. 101, 132501 (2008). 

[64] A. Gezerlis and J. Carlson, Phys. Rev. C81, 025803 
(2010). 



6 



[65] I. Vidana, A. Polls, and A. Ramos, Phys. Rev. C65, 
035804 (2002). 

[66] A. W. Steiner, J. M. Lattimer, and E. F. Brown, Astro- 

phys. J. 722, 33 (2010). 
[67] F. J. Fattoyev, C. J. Horowitz, J. Piekarewicz, and 

G. Shen, Phys. Rev. C82, 055803 (2010). 
[68] H. Mueller and B. D. Serot, Nucl. Phys. A606, 508 

(1996). 

[69] B. G. Todd and J. Piekarewicz, Phys. Rev. C67, 044317 



(2003). 

[70] F. Ozel, G. Baym, and T. Giiver, Phys. Rev. D82, 
101301 (2010). 

[71] F. Ozel, D. Psahis, R. Narayan, and A. S. ViUarreal, 

Astrophys. J. 757, 55 (2012). 
[72] J. Abadie et al, (LIGO Scientific Collaboration, Virgo 

Collaboration), Class. Quant. Grav. 27, 173001 (2010). 



